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Abstract 

We present a surface stacking fault (SSF) energy approach to predicting defect nucleation from the surfaces of surface-dominated 
nanostructure such as FCC metal nanowires. The approach is validated through a comparison between the SSF energy calculation 
and low-temperature classical molecular dynamics simulations of copper nanowires with different axial and transverse surface 
orientations, and cross sectional geometries. We focus on the effects of the geometric cross section by studying the transition from 
slip to twinning previously predicted in moving from a square to rectangular cross section for (100)/{100} nanowires, and also for 
moving from a rhombic to truncated rhombic cross sectional geometry for (110) nanowires. We show that the SSF energy criteria 
successfully predicts the initial mode of surface-nucleated plasticity at low temperature, while also discussing the effects of strain 
and temperature on the applicability of the criterion. 
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1. Introduction 

One of the fundamental mechanical properties of a material 
is its yield strength, or the maximum amount of load a material 
can take before sustaining irreversible plastic deformation, and 
eventually failure or fracture. In metals, because the unit plastic 
deformation mechanism is the dislocation, researchers have fo- 
cused on developing models to predict when a metal will yield 
via dislocation nucleation, and therefore what the yield strength 
of the metal will be. 

Seminal work towards this end was performed by Vitek 
Ollll], wno formulated a curve of potential energy versus dis- 
placement known as the generalized stacking fault (GSF) curve. 
The curve typically takes a nearly sinusoidal shape j3], in which 
the energy barrier that must be crossed in order for defect nu- 
cleation to occur is called the unstable stacking fault energy y us 
(see Fig. [TJ. Upon crossing the y us barrier, an equilibrium lat- 
tice spacing that does not correspond to the bulk equilibrium 
structure can be found, which corresponds to an intrinsic stack- 
ing fault (ISF). The resulting excess energy as compared to the 
bulk equilibrium for this configuration is called the intrinsic 
stacking fault energy ji s f. In the case of face-centered cubic 
(FCC) metals, the displacement axis of the GSF curve is along 
the (112) direction as this is the slip direction for partial dislo- 
cations in an FCC crystal. 
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The ideas underlying the GSF have been utilized for various 
predictions regarding nanoscale plasticity over the years. One 
important example is the work of Rice |4j], who formulated a 
model of dislocation nucleation from a crack tip. Specifically, 
Rice considered the competition between crack tip yielding via 
dislocation nucleation as compared to brittle fracture by com- 
paring the energy needed to create new fracture surfaces to the 
energy barrier to be overcome (y us ) in order for ductile fracture 
to occur. Rice's ideas were extended by Tadmor and Bernstein 
and Bernstein and Tadmor [6] to consider the competition 
between dislocation nucleation and twinning at a crack tip. 

There has been considerable success, as demonstrated above, 
in applying the GSF to predicting, at an atomistic level, defect 
nucleation in bulk materials. However, the validity of the GSF 
for predicting defect nucleation from the surfaces of nanostruc- 
tures such as nanowires has not been demonstrated to-date. 

Nanowires are an important class of low-dimensional nanos- 
tructures that have been extensively studied in the past ten or 
so years 10]. Because they are one-dimensional in nature, and 
because they are characterized by a large surface area to vol- 
ume ratio, their transverse surfaces have been found to play a 
critical role in controlling their operant inelastic deformation 
mechanisms. The deformation mechanisms of metal nanowires 
have been studied in detail recently, primarily through classical 
molecular dynamics (MD) simulations. While more extensive 
reviews on the deformation and plasticity of metal nanowires 
have recently appeared iS @], we highlight here those works 
that are particularly relevant to the present work, i.e. those 
which demonstrate the unexpected surface-mediated inelastic 
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deformation mechanisms in nanowires that motivate the present 
work. 

One important finding was that the transverse surfaces have 
a direct, first order effect on the operant inelastic deforma- 
tion mechanism ifioll of metal nanowires. This work lioll also 
demonstrated that traditional notions used to predict deforma- 
tion bulk mechanisms such as the Schmid factor may not be 
able to capture the surface-dominated mechanisms observed in 
nanowires. Another important finding was that the geometry of 
the nanowire cross section also has a direct effect on the oper- 
ant inelastic deformation mechanism of nanowires ll l[ll2lll3ll . 
Ji and Park 111 ll Il2ll elucidated this by considering the defor- 
mation of (100) copper nanowires with cross sectional geome- 
tries varying from square to rectangular. They found that while 
square nanowires yielded by nucleation and propagation of full 
and partial dislocation slip, rectangular nanowires yielded by 
nucleation and propagation of twins. A similar conclusion was 
found by Leach et al. II 1 311 . who studied the effects of truncat- 
ing a rhombic cross section on the deformation mechanisms 
of silver nanowires. While fully rhombic cross section (110) 
nanowires were observed to reorient in tension via nucleation 
and propagation of twins to a (100) orientation lfl4l [T^l . trun- 
cated rhombic nanowires were found to fail via nucleation and 
propagation of partial dislocations, while pentagon cross sec- 
tion nanowires were also found to fail via nucleation and prop- 
agation of partial dislocations. 

These MD studies collectively demonstrate that a predic- 
tive model of surface-nucleated incipient plasticity should ac- 
count for loading mode (i.e. tension or compression), trans- 
verse surface orientation, nanowire axial orientation and the ge- 
ometry of the nanowire cross section. It is worth noting that 
while the experimental literature on deformation mechanisms 
in metal nanowires is quite sparse compared to the fully popu- 
lated literature on MD simulations, evidence of the effects dis- 
cussed above has emerged. For example, the brittle nature of 
deformation and fracture of pentagonal silver nanowires was 
recently demonstrated by Zhu et al. Jlal . Similarly, rhom- 
bic cross section (110) gold nanowires were observed to un- 
dergo long-ranged twin migration and propagation, enabling 
a tensile-stress-induced reorientation to square cross section 
(100) nanowires lfl7ll . In contrast, pristine hexagonal cross sec- 
tion copper nanowires were observed to undergo brittle fracture 
by Richter et al. JlUl . which suggests that the cross sectional 
geometry of the nanowires is in fact a critical factor controlling 
the inelastic deformation mechanism. 

Outside of classical MD simulations, there have also been 
atomistic approaches based on classical stability theory, in 
which defect nucleation at the surfaces of nanostructures oc- 
curs when either the determinant of the Hessian becomes neg- 
ative [19], or when the change in potential energy from one 
configuration to another decreases 12011 . However, both of these 
approaches were for defect nucleation from infinitely large sur- 
faces attached to semi-infinite bodies, and therefore may not 
hold predictive power for defect nucleation from the surfaces 
of one-dimensional nanowires. Other researchers have also uti- 
lized energy landscape exploration techniques to investigate the 
nucleation barriers that must be overcome for specific plas- 



tic deformation mechanisms to occur These ap- 

proaches, while insightful, have also not yielded a generic and 
predictive capability that correlates the observed incipient plas- 
tic deformation mechanism to the nanowire size, cross sectional 
geometry, axial and surface orientation. 

Therefore, we present in this work a generalized stacking 
fault criteria for surfaces. The work is based upon the classical 
GSF approach of Vitek 01 though applied to finite domain 
nanostructures such as nanowires. Ideas very similar to the ones 
presented in this manuscript were discussed by Leach l24ll in a 
Ph.D. thesis, though no work along these lines was published. 
We demonstrate that because this approach accounts for all fea- 
tures discussed above, i.e. loading mode, transverse surface 
orientation, axial orientation and cross sectional geometry, it is 
applicable to a wide range of metal nanowire geometries and 
loading modes. We illustrate the proposed approach as com- 
pared to companion benchmark MD simulations. 

2. Methodology 

The interaction between copper atoms is described using the 
embedded atom method (EAM) 112.511 - using the potential devel- 
oped by Mishin et al. 112611 . In the EAM, the total energy of the 
system is expressed in terms of the summation over each atom 



E = Y j F i (p i ) + Y J <f> i j(r ij ) > 



(1) 



where N is the total number of atoms. Atom ; is embedded in its 
position by the first 'embedding' term with respect to the elec- 
tron density p,- at this position. The second term provides the 
repulsive force between atom ; and its neighbors, which helps 
to stabilize the whole structure. The EAM potential of Mishin 
et al. i26ll was chosen due to its accuracy in capturing the stack- 
ing fault and surface energies of copper. 

All MD simulations were performed using the publicly avail- 
able simulation code LAMMPS JI^Hl, while the OVITO 
package was used for visualization I29ll. In the MD simula- 
tion, the system is first thermalized at a constant temperature by 
the Nose-Hoover ll30L 31] heat bath for 100 ps. Free boundary 



conditions were applied along the two lateral directions of the 
nanowire by introducing a large vacancy space in the simula- 
tion box along these two directions. Along the axial direction, 
one end of the nanowire was fixed while the other was free. 



Tensile strain was then applied, similar to previous works [32], 
by adding a ramp velocity that scaled linearly from zero at the 
fixed end to a maximum value at the free end. The added ve- 
locity profile results in a strain rate of e = 10 9 s -1 , which is a 
typical value in MD simulations, and which does not preclude 
the generation of twins, full or part ial dislocations, as shown in 
previous works OS [II El Bill El- 

For bulk FCC copper, the GSF energy curve is calculated us- 
ing standard procedures 01 . Specifically, the simulation block 
is rectangular with its three orthogonal coordinate directions 
along the (111), (110), and (112) lattice directions. Periodic 
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Figure 1 : (Color online) The GSF energy curve for bulk FCC copper for slip on 
a j 1 1 1 ) plane along the (112) direction. The experimental value for the intrinsic 
stacking fault energy is shown as a dot at point A. The three extreme points are 
labeled according to the three corresponding configurations shown in the top 
inset. The unstable stacking fault energy (y us ) and the intrinsic stacking fault 
energy (y; s f) is displayed. 



boundary conditions are applied in the (110) and (112) direc- 
tions, while free boundary conditions are applied in the (111) 
direction. The structure is divided into two halves by a {111} 
plane. These two halves are sheared with respect to each other 
along the (112) direction using small displacement increments. 
The total potential energy (E tot ) of the structure is recorded at 
each shearing displacement, while the atoms are allowed to 
relax in the (111) direction during the shearing process. The 
GSF energy is then calculated from (E tot - £^ t )/A(ni}, where 
(£ to t _ £ t °ot) i s me change in the total potential energy due to the 
relative displacement between two halves and A{ui) is the area 
of the shearing {111} plane. 

Fig.[TJshows the GSF energy curve for bulk copper. The top 
inset in the figure demonstrates the evolution of the atom con- 
figuration at the shearing interface during the relative displace- 
ment of the two halves of the crystal. An intrinsic stacking fault 
(ISF) occurs at the partial Burgers vector b p with b p — a/ y6, 
while an unstable stacking fault (USF) occurs at the displace- 
ment of b p /2. As expected, good agreement in the ISF energy 
jisf is obtained between the Mishin et al. 12611 potential (44.5 
mJ/m 2 ) and experiments (45.0 mJ/m 2 ) l33ll . 

The surface GSF energy curve for unstrained nanowires can 
be calculated in following two steps. First, the nanowire is 
relaxed to a minimum energy configuration. During this pro- 
cess, surface atoms contract inwards towards the bulk atoms 
to maximize their electron density. Second, the surface GSF 
energy curve for the twinning fault of the nanowire can then 
be obtained through the process shown in Fig. |2 where the 
(100)/{100) CuNW in (a) is sheared rigidly along a { 1 1 1 ) plane 
in the (112) direction. Unlike in the bulk case, the two blocks 
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Figure 2: (Color online) The formation of a twinning fault in a (100)/{100) 
CuNW, where the (111) plane slips in the (112) direction. 
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of atoms comprising the nanowire are not allowed to relax in 
the nanowire case. This is primarily due to the difficulty in 
constraining the nanowire relaxation to occur only in the (111) 
direction, as in the bulk case. However, as we demonstrate be- 
low, this small approximation does not degrade the accuracy or 
utility of the surface GSF curves. As the strain increases, the 
configuration of the unstable stacking fault is reached in (b), 
followed by the formation of an ISF in (c). After the first ISF is 
formed, the deformation continues to shear a neighboring {111} 
plane in the same (112) direction, resulting in the formation of 
the second unstable twinning fault seen in (d). Finally, the for- 
mation process finishes with the twinning fault in (e). Similar 
to the bulk case, the surface GSF energy for nanowire here is 
then calculated from (E tot - £^ t )/A(in}. 

We also calculate the surface GSF energy curves for axially 
strained nanowires, for reasons that we explain in detail in the 
discussion of the numerical examples. For these surface GSF 
curves, the nanowire is strained axially in increments of 0. 1 %, 
where the atoms are allowed to relax after the two ends of the 
nanowire are fixed at the new length. From these energy min- 
imizing positions, the shearing process described above and il- 
lustrated in Fig. |2]is performed to extract the surface GSF en- 
ergy curve for the prescribed axial strain. 

Fig.[3]demonstrates the formation process of a slip fault. The 
production of the first ISF from (a) to (c) is the same as the 
formation process of the twinning fault shown in Fig.|2j so the 
surface GSF energy curve are the same for twinning and slip 
faults up to the occurrence of the first ISF in Fig. 0c). The 
slip formation continues by generating a second ISF, where the 
nanowire is sheared along another {111} plane along its corre- 
sponding (112) direction. The slip fault is formed in (e). The 
surface GSF energy for the slip fault is also calculated from the 
change of the total potential energy in the structure. 

The surface GSF energy curves for twinning and slip defor- 
mation pathways are the same up to the emergence of the first 
ISF. Afterwards, these two GSF energy curves start to dispart. 
At the second unstable stacking fault, if the twinning surface 
GSF energy is lower than that of the slip deformation, then the 
twinning deformation will be the energetically favored pathway 
during the inelastic deformation of this system. Otherwise, slip 
will be the energetically favored pathway. In this way, the sur- 
face GSF energy curve can be used as a criteria to predict the 
energetically favorable plasticity mechanism during the defor- 
mation of the nanowire. 




Unstable 
Stacking Fault 



(c) 
Intrisic 

Stacking Fault 



Unstable 
Slip Fault 




Figure 3: (Color online) The formation of a slip fault in a < 100>/{ 100) CuNW. 



3. Results 

3.1. Square to Rectangular NWs 

The rationale for the simulations in this section are to inves- 
tigate the ability of the proposed surface GSF criteria to cap- 
ture the transition from slip to twinning-dominated deforma- 



tion that was observed in earlier MD simulations 111 111 when 
the cross section of (100)/{100) copper nanowires (CuNWs) 
changed from square to rectangular. We perform such simula- 
tions, while calculating the surface GSF energy curves for both 
twinning and slip deformation pathways. These surface GSF 
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Figure 4: (Color online) The surface GSF energy curve for three CuNWs with 
size (l x ,l y ,l z )=(40, 4, 4), (40, 4, 8), and (40, 4, 12). Curves for twinning and 
slip deformation are displayed by solid and dashed lines, respectively. Lines 
are colored according to the color of the cross section in the inset in (a). The 
surface GSF energy for structures with strains 0.0, 0.02, and 0.05 are shown in 
panels (a), (b), and (c). 



energy curves are then used to predict for the favored defor- 
mation pathway in the CuNWs. Concurrent low-temperature 
(0.1K) MD simulations are performed to monitor the deforma- 
tion pathway during the inelastic deformation of the CuNWs 
under tension. We will validate the surface GSF criteria by 
comparing the deformation pathway predicted by the surface 
GSF energy calculations to that observed in the MD simula- 
tions. 

The three nanowires studied had size (l x , l y , y =(40, 4, 4), 
(40, 4, 8), and (40, 4, 12) cubic lattice units (CLU), where the 
CLU=a = 3.615 A is the lattice constant of FCC copper. All 
three CuNWs had the same axial orientation along the (100) 
direction, with four {100} side surfaces. Fig. |4](a) shows the 
surface GSF energy curves of the undeformed CuNWs. For the 
two CuNWs with l z = 4 and 8, the deformation energy required 
for slip is lower than that for the twinning. For L = 12, the 
surface GSF energy curves for the twinning and slip pathways 
are almost indistinguishable, which implies that either twinning 
or slip could occur. 
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Figure 5: (Color online) The configuration for the CuNW (40, 4, 4) from MD 
simulation at 0.1 K with strain e = 0.17. Atoms are colored by potential energy 
in the top image and by the centrosymmetry parameter in the bottom image. 



We also show the surface GSF energy curves for all the three 
CuNWs with strain e = 0.02 and 0.05 in Figs. g](b) and (c), 
where we emphasize that these strains are below the yield strain 
of all CuNWs. By performing the surface GSF calculation on 
the tensile-strained nanowires, it is clear that the strain has no 
effect on the square CuNW (40,4,4), where slip remains the 
preferred deformation mechanism. 

In contrast, the difference in the energy barrier for twinning 
and slip in the CuNW (40, 4, 8) becomes very small at e = 
0.05, with the surface GSF energy for slip being slightly lower 
than twinning. We have also calculated the surface GSF energy 
curve for CuNW (40, 4, 8) with even larger tensile strains, and 
observed similar trends, i.e that the energy for slip is always 
very close to but slightly lower that for twinning. However, as 
the cross sectional geometry becomes even more rectangular 
as for the CuNW (40,4,12) case, it is seen that the twinning 
pathway energy barrier becomes progressively smaller as strain 
increases than that for slip. These results are in agreement with 
the MD simulation results of Ji and Park lUlll . who predicted 
the slip to twinning transition as the nanowire cross sectional 
geometry transitioned from square to rectangular. 

We verify the predictions of the surface GSF curves in Fig. [4] 
by performing MD simulations of the tensile deformation of 
CuNWs. Because of the strain-induced evolution in the pre- 
dicted deformation mechanism seen in Fig. [4] we use the sur- 
face GSF energy calculated from a strained CuNW in Fig. HJc) 
to predict the deformation pathway we expect to see in the 
MD simulations. From Fig. |4](c), the energetically preferred 
deformations predicted by the surface GSF energy are slip, 
slip/twinning, and twinning for the three CuNWs with l z — 4, 
8, and 12. For I- = 8, two deformation pathways can occur si- 
multaneously due to their close surface GSF energy curves, and 
thus our surface GSF calculation predicts a transition from slip 
to twinning deformation in CuNWs (l x , l y , l z ) with increasing 
width to height ratio L/l y . 

Fig.|5]shows the configuration of the CuNW (40, 4, 4) under 
strain e = 0.17 at 0.1 K. The top image is colored according to 
the potential energy per atom, while the bottom image is col- 
ored by the centrosymmetry parameter. Only atoms with cen- 
trosymmetry parameters within the color bar are shown in the 
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bottom image. No twinning deformation is observed in the de- 
formed structure, while slip deformations can be observed. This 
MD simulation result is therefore consistent with the prediction 
by the surface GSF energy in Fig.|4](c). 

Fig.[6]shows the evolution of the configuration for the CuNW 
(40, 4, 8) from MD simulation at 0. 1 K. With increasing ten- 
sion, the CuNW starts to yield at e = 0.104, where the initial 
yield mechanism is partial dislocation slip from the top and bot- 
tom surfaces that meet at the center of the nanowire to form a 
rhombic dislocation structure as illustrated by the bottom im- 
age for e = 0.104. Under slightly higher tension, e = 0.107, the 
rhombic deformation shapes break open and propagate along 
the axial direction of the nanowire. At e = 0.15, a transition 
from slip to twinning is observed. These MD simulation re- 
sults are consistent with the surface GSF energy prediction in 
Fig.ffl(c). 

Fig.|7]shows the deformation of the CuNW (40, 4, 12) under 
tensile loading at 0.1 K. Twinning is found to be the sole de- 
formation mechanism throughout the MD simulation. Specifi- 
cally, there are several observable twin boundaries in separate 
regions of the nanowire at e = 0.11. With increasing strain 
e = 0.2, these twin boundaries separate to form larger twinned 
regions in the nanowire. The MD simulations again verify the 
prediction from the surface GSF energy in Fig.|4](c). 

Overall, we have shown in this section that the surface GSF 
criteria can correctly predict the previously observed transition 
from slip to twinning in CuNWs as the cross sectional geometry 
evolves from square to rectangular, particularly if strain effects 
on the surface GSF curve are accounted for. 

3.2. Rhombic NWs 

The rationale for the simulations in this section are to inves- 
tigate the ability of the proposed surface GSF criteria to capture 
the transition from slip to twinning-dominated deformation that 
was observed in earlier MD simulations H 1 3fl when the cross 
section of(110)/{lll) silver nanowires changed from rhombic 
to truncated rhombic. We investigate this phenomenon in the 
present work for copper nanowires (CuNWs). The truncations 
considered are shown in Fig. [8] where the cross section is trun- 
cated with truncation level n = 0, 1, 2, 3, and 4. For n — 0, 
all side surfaces are {111} planes, with the length of this set of 
CuNWs being 10 nm. The CuNW is truncated at the two lower 
apex, resulting in the exposure of two {100} side surfaces. With 
increasing truncation level n, the area of the exposed {100} sur- 
faces increases. In the following, we will first present the sur- 
face GSF energy calculation for this set of CuNWs, followed 
by the MD simulation results. 

The surface GSF energy curves for all five truncated rhombic 
CuNWs are shown in Fig. [9] (a) for structures without strain. 
The curves for twinning are displayed by solid lines while the 
curves for slip are displayed by dashed lines. All curves are 
colored according to the corresponding cross section as shown 
in the inset of Fig. [9] (a). With increasing truncation levels n 
from to 4, there is a transition in the relationship between 
the deformation energy of twinning and slip pathways. The 
transition happens at n = 2, where the energy barrier difference 
between twinning and slip pathways is almost zero. 
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Figure 6: (Color online) The evolution of the configuration for CuNW (40, 
4, 8) at 0. 1 K with increasing tensile strain. The top image in each figure is 
colored according to potential energy while the bottom image is colored by the 
centrosymmetry parameter. At e = 0.104, partial dislocations from the top and 
bottom surfaces form two rhombic shapes. These rhombic shapes propagate 
along the CuNW under increasing tensile strain before e < 0.107 and turn into 
twinning deformation at e = 0.15. 
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Figure 7: (Color online) Tensile deformation of the CuNW (40, 4, 12) at 0.1 K. 
Twinning is the only observed plastic deformation mechanism. 




Figure 8: (Color online) Cross section of the rhombic CuNW with truncation 
level n = 0, 1, 2, 3, and 4. The side view in the last image displays the (110) 
growth direction of the CuNW. For perfect rhombic CuNW, all side surfaces 
are j 1 1 1 1 planes. Two j 100) lateral side surfaces will be exposed by truncation. 
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Figure 9: (Color online) The surface GSF energy curves for truncated rhombic 
CuNWs. Curves are displayed by solid lines for twinning and by dashed lines 
for slip deformation. Lines are colored according to the color of the cross sec- 
tion in the inset in (a). Panel (b) shows the surface GSF energy curve around 
the second maximum region for nanowires with strain e = 0.05. 
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Figure 10: (Color online) The evolution of the rhombic < 1 10>/{1 1 1 } CuNW 
with truncation level n = from MD simulation at 0. 1 K. Atoms are colored 
by the centrosymmetry parameter. Reorientation of the initially rhombic (1 10) 
CuNW to a square cross section (100)/{100) CuNW via long-range twin bound- 
ary migration is observed. 



As discussed above for the square to rectangular cross sec- 
tion nanowires, accounting for the effects of strain is impor- 
tant to accurately apply the surface GSF energy criteria. We 
thus show in Fig. [9] (b) the surface GSF energy for these five 
CuNWs under tensile strain e — 0.05. For clarity, only curves 
around the second extreme values in Fig. [3 a) are shown. We 
will use Fig.|9](b) to predict the preferred deformation pathway 
for the truncated rhombic CuNWs. From these surface GSF en- 
ergy curves, the deformation energy for twinning is lower than 
that of slip in rhombic CuNWs with truncated level < « < 2. 
As a result, the twinning deformation is predicted to be the pre- 
ferred pathway of these three CuNWs. The deformation energy 
for twinning is higher than that of the slip in rhombic CuNWs 
with truncated level n > 3. It is thus predicted that the slip 
deformation is the preferred pathway in these two (n = 3,4) 
CuNWs. We will perform MD simulations to show the transi- 
tion from twinning to slip in the deformation pathway for this 
set of CuNWs with increasing n from to 4. Similar MD sim- 
ulations can be found in the work by Leach et.al lfl3ll . 

Fig. [TOl shows the tensile deformation at 0.1 K in the MD 
simulation for the rhombic CuNW with truncation level n — 0. 
The only observed deformation mechanism is twinning, which 
is generated at lower strain and then propagates along the 
nanowire length under further tensile loading. The long-range 
propagation of the twin boundaries reorients the side surfaces 
from {111} to {100}, while the axial orientation reorients from 
(110) to (100) with a corresponding change in cross sectional 
geometry from rhombic to square. These simulation results are 
consistent with earlier MD simulations focused on nanowire 
shape memory and pseudoelasticity Il4[ll5[ 13411 . and more re- 
cently the work of Seo et al. [vh . who gave the first experimen- 
tal demonstration of such long-ranged twin boundary propaga- 




Figure 1 1 : (Color online) The evolution of the rhombic CuNW with truncation 
level n = 3 from MD simulation at 0. 1 K. Atoms are colored according to the 
centrosymmetry parameter. At the initial deformation stage, a slip deformation 
is generated at e = 0.08, where the dislocated atoms are removed only in this 
particular image to clearly show the slip deformation. Under further tension, 
the deformation pathway switches to twinning. The twinning deformation then 
propagates along the entire length CuNW. 



tion and the subsequent reorientation. 

After the propagation of the twinning boundaries along the 
nanowire length, the reorientated (100}/{100) nanowire under- 
goes yielding via slip at e = 0.5. We similarly find through the 
MD simulations that twinning is the only deformation mech- 
anism observed during the tensile deformation of the rhombic 
CuNWs with truncated levels n — 1 and 2. These MD simula- 
tion results are in good agreement with the surface GSF energy 
shown in Fig.[9](b). 

Fig.QT|shows the MD simulation for the tensile deformation 
of the rhombic CuNW with truncation level n = 3 at 0.1 K. 
In the second image, copper atoms within the deformed region 
are removed, so that a slip deformation can be clearly observed 
near the right boundary at strain e = 0.08. With increasing 
tensile strain, this slip deformation turns into a twinning de- 
formation as shown in the third image at e = 0.1. The twin- 
ning deformation is then propagated along the nanowire and 
the CuNW is reorientated from a rhombic (110) CuNW into a 
square (100) nanowire with two {100) surfaces and two {110) 
surfaces. We note that the reorientation of the side surfaces in 
going from the rhombic (110) to square (100) nanowire is as 
follows: {111) — > {100) and {100) — -> (110). An important 
point to make here is that the MD simulations demonstrate that 
while the surface GSF energy calculation accurately predicts 
the initial deformation mechanism in the nanowires, it does not 
yield accurate predictions of subsequent plastic events at larger 
strains. This is not surprising considering that the surface GSF 
calculation is done for defect-free, strained structures without 
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Figure 12: (Color online) The evolution of the CuNW with truncation level 
n = 4 from MD simulation at 0. 1 K. Atoms are colored by the potential energy 
in the top image and by the centrosymmetry parameter in the bottom image. 
Slip is the only deformation mechanism observed at this temperature. 
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Figure 13: (Color online) The evolution of the CuNW with truncation level 
n = 4 from MD simulation at 10 K. Atoms are colored by the centrosymme- 
try parameter. Twinning is the only observed deformation mechanism at this 
temperature. 



any pre-existing stacking faults or twins. 

Fig. Q~2] shows the slip-dominated deformation of the rhom- 
bic CuNW with truncation level n — 4 from MD simulation at 
0. 1 K, where the top image is colored according to the potential 
energy per atom while the bottom image is colored by the cen- 
trosymmetry parameter. The slip deformation and the absence 
of any twinning is clearly observed in the bottom image. The 
MD simulation is consistent with the surface GSF energy cal- 
culation, which predicts the slip deformation to be the preferred 
pathway in this CuNW. 

To study temperature effects on the deformation pathway, we 
show in Fig.QjJthe structure evolution from MD simulation for 
the truncated rhombic CuNW with n — 4 at 10 K. It is clear that 
twinning is the only deformation mechanism that is observed, 
where the twinning leads to the same reorientation previously 
observed, for example, in Fig. [10] and in general for all CuNWs 
with truncation levels n < 2 at low temperature 0. 1 K. 

The MD simulation at 10 K thus illustrates that the surface 
GSF energy calculation cannot predict the deformation path- 
way at higher temperature, which can be understood from an 



energetic perspective. Fig. [9] (b) shows that the energy bar- 
rier between twinning and slip is about AE = 14.2 mJ/m 2 at 
the second unstable stacking fault. Most of this energy comes 
from those copper atoms sitting in the two {111} planes that 
are sheared along the (112) direction to calculate the surface 
GSF energy. The energy barrier can further be distributed onto 
each copper atom in these two {111} planes by multiplying the 
area per atom in this {111} plane: i.e ^gsf = A£ x Ao/2, with 
Aq = V3a 2 /4 = 5.7 A 2 , where the factor of 2 is needed to 
divide the energy barrier onto one of the two shearing {111} 
surfaces. Some simple algebra then gives £gsf = 2.5 meV. 

Besides £gsf, there is another important energetic contribu- 
tion in this system at finite temperature, i.e the thermal vibra- 
tion energy. According to the equipartition theorem, the vibra- 
tional energy of each degree of freedom at temperature T is 
ksT/2, where ks is the Boltzmann constant. The thermal en- 
ergy for each atom is Ej = 3^772, which is 1.3 meV at 10 K. 
We can see that the two energy scales £gsf and Ej are on the 
same order. The thermal vibration reflects a random movement 
of each copper atom, while the GSF energy curve predicts a 
deterministic movement. Since the energy scales of these two 
movements are on the same order, the random movement from 
thermal vibrations can counteract the deterministic motion im- 
posed in calculating the surface GSF energy, and is the reason 
why the surface GSF prediction fails at higher (i.e. non-zero) 
temperature. In contrast, at the low temperature 0.1K that we 
have used for the majority of our MD simulations, the thermal 
vibrational energy Ej is about two orders smaller than the sur- 
face GSF energy barrier £gsf, so the thermal-induced random 
movement can be safely ignored, resulting in accurate predic- 
tions from the surface GSF energy. 

Before concluding, we provide an energetic explanation for 
the twinning to slip deformation transition in the truncated 
rhombic CuNWs. During the deformation of the truncated 
rhombic CuNWs, the four {111} surfaces are reoriented to {100} 
surfaces, while the two {100} surfaces are reoriented to {110} 
surfaces. Furthermore, it is important to note that in both 
cases, the reoriented surface has a higher surface energy than 
the initial surface orientation. Because of this, we calculate 
the surface energies for the {111}, {100} and {110} surfaces to 
be 1242.3, 1352.6, and 1480.5 mJ/m 2 , which are in agreement 
with the values published by Mishin et al. l26ll . Based on these 
values, it can be seen that there is an energetic penalty of 109.3 
mJ/m 2 to reorient a {111} surface into a {100} surface. Simi- 
larly, there is a 127.9 mJ/m 2 energetic penalty for reorienting a 
{100} surface into a {1 10} surface. In both cases, these energetic 
differences represent energy barriers that must be overcome by 
the twinning deformation. As the truncation level increases, a 
greater proportion of {100} surface area is exposed with a sim- 
ilar reduction in {111} surface area. This is important because, 
as discussed above, the energetic cost of reorienting a {111} 
surface to {100} is lower than the cost of reorienting a {100} 
surface to {110}. This is why the energetic barrier to twinning 
deformation increases as the truncation of the rhombic CuNWs 
increases. 

In contrast, when slip is the operant deformation mode, 
no lattice reorientation is observed in the truncated rhombic 
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CuNWs. In this case, the energy that is added to the system 
results from the imperfect stacking arrangement in the NW in 
conjunction with the surface step that results at the free surface 
due to the stacking fault. The energetic barrier for this process 
also changes with increasing truncation, but not as dramatically 
as it does during the twinning-dominated deformation mode. 
As a result, there is a predicted twinning to slip transition in the 
deformation mechanism as discussed above. 

4. Conclusions 

In conclusion, we have proposed a criteria based on the 
surface stacking fault (SSF) energy to predict the incipient 
mode of plasticity in surface-dominated nanostructures such 
as nanowires. The validity of the idea was demonstrated by 
comparison with classical molecular dynamics simulations of 
nanowires that were previously predicted to exhibit a deforma- 
tion mechanism transition from slip to twinning based on vari- 
ations of the nanowire cross sectional geometry. The SSF cri- 
teria was shown to be most accurate at low temperatures, while 
strain effects were demonstrated to potentially alter the defor- 
mation mechanism from slip to twinning for certain nanowire 
geometries. 
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